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The resistance to fracture of reversible biopolymer hydrogels is an important 
control factor of the cutting/slicing and eating characteristics of food gelsi. It is 
also critical for their utilization in tissue engineering, for which mechanical pro- 
tection of encapsulated components is needed2.j2i. Its dependence on loading rate^ 
and, recently, on the density and strength of cross-links^ has been investigated. 
But no attention was paid so far to solvent nor to environment effects. Here we 
report a systematic study of crack dynamics in gels of gelatin in water/glycerol 
mixtures. We show on this model system that: (i) increasing solvent viscosity 
slows down cracks; (ii) soaking with solvent increases markedly gel fragility; (iii) 
tuning the viscosity of the (miscible) environmental liquid affects crack propaga- 
tion via diffusive invasion of the crack tip vicinity. The results point toward the 
fact that fracture occurs by viscoplastic chain pull-out. This mechanism, as well 
as the related phenomenology, should be common to all reversibly cross-linked 
(physical) gels. 

Gelatin gels are constituted of denatured (coil) collagen chains, held together by cross- 
links made of segments of three-stranded helices stabilized by hydrogen bonds^. This net- 
work, swollen by the aqueous solvent, which controls its (undrained) bulk modulus, is re- 
sponsible for the finite shear modulus /i, of order a few kPa. Hence, hydrogels can be 
considered incompressible. One estimates average mesh sizes ^ ~ (A;T//i)i/3 of order 10 nm, 
i.e. coil segments involving a few 100 units (residues)^. Moreover, in the presence of pressure 
gradients, the solvent diffuses through the network. This poroelastic behaviour'''^ controls 
e.g. slow solvent draining in or out of the gel under applied stresses. 

They are thermoreversible, i.e., in contrast with chemical, covalently cross-linked gels, 
their network "melts" close above room temperature. This behavior, assignable to their 
small cross-link binding energy, leads to the well studied^ slow aging (strengthening) of /i, 
and to their noticeable creep under moderate stresses^. When stretched at constant strain 
rate, gelatin gels ultimately fail at a strain ~ 1 which, though rather poorly reproducible, is 
clearly rate-dependenl4. In order to get insight into the nature of the dissipative processes at 
play, one needs to investigate the propagation of cracks independently from their (stochastic) 
nucleationi^. Here we study the fracture energy G{V) needed to propagate a crack at 
constant velocity V in notched long thin plates (see Fig. of gels differing by the glycerol 
content of their aqueous solvent. 
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As seen on Figure I21^a), Q is very strongly velocity and solvent-dependent. Our gels 
are velocity toughening: at fixed (p, Q grows quasi-linearly over the whole investigated 
V^-range (0.1-30 mm.s^^). Linear extrapolation down to V = yields an evaluated quasi- 
static toughness Qq. Within experimental accuracy, Qq is of order 2.5 J.m~^, i.e. about 20 
times larger than a gel-air surface energy, and 0-independent. In contrast, QiV) becomes 
noticeably steeper as increases. Plotting it versus rjV (Fig. Efb)) captures most of this 
dependence. Note that the ratio Q/rjV is a huge number of order 10^. 

This points toward the critical role of network/solvent relative motion. In particular, the 
impossibility for not very thin, quasi incompressible plates to accomodate fully the high local 
strain gradients developing in the crack tip region results in high negative pressures. So, 
most likely, solvent is partly drained out of this region into the bulk, leading for the chains 
thus exposed to air, to a solvation energy cost. We investigate this issue with the help of 
experiments in which a drop of the gel solvent is introduced into the already moving crack 
opening. Such tip wetting induces, at fixed sample stretching, a positive, ^-independent 
velocity jump (Fig. IH)). Equivalently, wetting decreases G{V) by a constant A^o- For = 0, 
AQq ~ 2 J.m^^ is a substantial fraction of Qq. 

Can one make, on the basis of these results, a plausible guess about the nature of the 
fracture mechanism in reversible hydrogels? 

Clearly, one cannot invoke here the classical Lake-Thomas picture'''^, which successfully 
accounts for rubber toughness. Indeed, in these materials, fracture occurs via chain scission: 
polymer segments crossing the fracture plane are stretched taut until they store an elastic 
energy per monomer of order the covalent binding one, Uchain ~ a few eV. In thermoreversible 
gels, the corresponding force , ~ Uchain/ is more than two orders of magnitude larger than 
that, /* ~ Ucl/O'i which can be sustained by the i7-bond stabilized cross-links (CL). Ucl 
is the segmental unbinding energy introduced in the zipper-^ and reel-chaini^ models of gel 
elasticity. This leads us to postulate that, in the highly stressed "active crack tip zone", 
cross-links yield, up till the stretched chains are pulled out of the matrix. The threshold 
stress (J* at the onset of CL-yield can be estimated as cr* ~ /*/^^ = Ucl/O'^'^- With 
a ~ 0.3 nm, Ucl ~ 10^^ eV, ^ ~ 10 nm, a* ~ 500 kPa, two orders of magnitude larger than 
a small strain shear modulus. 

When solvent can be pumped from a wetting drop, the plastic zone deforms under this 
constant stress till the opening 5c at the tip reaches the length £ of a fully stretched chain. 
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This Dugdale-like picture^ yields for the quasi-static fracture energy : 
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(1) 



from which we get £ ~ 1.2 /xm. With an average mass M^es = 80 Da per residue, this yields 
for the gelatin molar weight the reasonable estimate: Mres^/ci ~ 320 kDa. 



Determining the value and precise space variation of a = d5/dx in the active zone would 
demand solving for the whole stress field. Given the large value of a* / fi, this would raise such 
intricate, still unsolved issues as strain-induced helix/coil transition^, elastic crack blunting 
and strain hardening at large deformation^^. Failing anything better, we take a to be an 
unknown geometrical factor. 

We then write for the viscoplastic stress a = a* + ^"^^(V^). The viscous stress a^s can 
be evaluated as resulting from hydrodynamic friction of chains of contour length i pulled 
at velocity 6 out of "tubes" formed by the embedding network. A natural candidate for 
the radius of these tubes is the effective pore size ^hydr extracted from light scattering 
experiments^^ (see Methods). So, ayis ~ ctV^^/^hydr^ 



So, this schematic model does account for the linear g{V) variation. Moreover, with 
^hydr = 2.7 nm, it predicts a reduced slope dQ /dirjV) ~ 2.5 x lO^a. When compared with 
experimental values (~ 10^) this suggests a values of order unity, possibly associated with 
elastic blunting^. 

Beyond this, the remaining splay between the Q{riV) curves (Fig.|2fb)) appears positively 
correlated with elastic stiffness variations: as is the case for reversible alginate gelsi^,^ "the 
stiffer, the tougher" — a relationship currenly under more systematic investigation. 

In this picture, draining at a non-wetted tip results in a capillary energy cost par chain 
At/o..^^ ~ 1000 eV/chain, i.e. ~ 10 fcT per residue, a value which suggests that chains are 
extracted individually rather than as gel fibrils. The observed V^-independence of [Q'^'^y — 
gwet^ indicates that a, hence the geometry of the active zone, is largely unaffected by wetting. 

We can now take one further step. Our scenario suggests that we should be able to 
tune, at fixed grips, the crack velocity by using a wetting drop with a glycerol content 
(pdrop, hence a viscosity, different from that of the gel solvent. We expect that, for small 



Let us turn to the ^-dependence of Q. A finite V means a finite pull-out velocity 6 = aV. 





4 



enough V, the miscible wetting fluid will invade the whole active zone, bringing V to its 
value for a 0dr-op-gel. The faster the crack, the less efficient this diffusive mixing process: the 
hetero- wetted G{V) should gradually approach that for the homo- wetted 0-gel. The result 
for cracks in a = 30% gel wetted by pure water, shown on Figure 4, spectacularly confirms 
this expectation. 

Moreover, assuming that the cross-over range {Vi,V2) (Fig. 0} corresponds to diffusion 
lengths Deff/V decreasing from dact to ^, we estimate for the diffusion constant of glycerol 
in the stretched gel -De// ~ ~ 2 x 10^^'' m^/s, and for the size of the active zone in the 
gel matrix dact ~ iV2/Vi ~ 100 nm. 

We conclude from this work that due to their weak binding responsible for cross-link 
plasticity^i^ii^, thermoreversible gels fracture via chain pull-out. While the fracture threshold 
is controlled by the cross-link yield stress, i.e. by the network only, crack dynamics is 
ruled by chain/solvent friction. This opens promising perspectives toward solvent control 
of crack dynamics in these materials, since: (i) The larger bulk solvent viscosity is, the 
slower cracks under a given loading, (ii) Homowetting of a crack tip speeds it up and can 
set subcritical precracks into motion, (iii) Heterowetting by a miscible fluid with substantial 
viscosity contrast leads, via diffusive "rinsing" of the active zone, to drastic effects on the 
propagation of slow cracks. 

So, tip wetting appears as a method of local and fast control of crack dynamics in such 
materials. Conversely, our results call attention to the fact that characterization of gel 
fracture properties should be performed under realistic, e.g. physiological, environmental 
conditions. 

METHODS 

Gel sample preparation 

Gels are prepared by dissolving 5 wt% gelatin powder (type A from porcine skin, Sigma) 
in mixtures of glycerol (0 = 0, 20, 30, 60 wt%) in deionized water, under continuous stirring 
at 90° C, an unusually high temperature needed to get homogeneous pre-gel solutions at 
high (j)^. A control experiment performed with a pure water/gelatin sample prepared at 60° 
C resulted in differences on low strain moduli and QiV) slopes of, respectively, 1% and 7%, 
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compatible with scatters between samples prepared at 90° C, whence we conclude that our 
preparation method does not induce significant gelatin hydrolysis. 

The solution is then poured into a mold consisting of a rectangular frame and two plates 
covered with Mylar® films. On the longest sides of the frame, the curly half-part of an 
adhesive Velcro® tape improves the gel plate grip. The mold is set at 2 ± 0.5° C for 15 hrs, 
then clamped to the mechanical testing set-up and left at room temperature (19 ±1° C) for 
1 hr. The removable pieces of the mold are subsequently taken off, leaving the 300 x 30 x 10 
mm^ gel plate fixed to its grips. The Mylar® films are left in position in order to prevent 
solvent evaporation. They are peeled off just before performing an experiment. 

Hydrogel characterization 

Small strain shear moduli fi are computed from the force-elongation curves assuming 
incompressibility, plane stress deformation and neglecting finite-size effects. A high glycerol 
content significantly stiffens the gels , from /i = 3.5 kPa at = to = 5.2 kPa at = 60 
%, possibly due to 0-dependent solvent-chain interactions. Large deformation, non-linear 
curves up to stretching ratios A ~ 1.5 are integrated numerically to determine the elastic 
energy ^(A) stored in stretched plates. 

The collective diffusive mode of the polymer network in the solvent^ is charaterized by 
the diffusion coefficient Dcoii, measured by dynamic light scattering as described elsewherei^. 
Solvent viscosities rj range from 10~^ Pa.s at = to 1 1 x 10^^ Pa.s at = 60 %. Accordingly, 
Dcoii decreases as increases. One estimates an effective pore size as $,hydr = a/ Dcourj/ fi. It 
is found independent of within experimental accuracy : ^hydr = 2.7 ± 0.2 nm. 

Fracture experiments 

Before stretching, a knife-cut notch is made at one edge of the plate. The grips are then 
pulled apart for 1 second by an amount Ah = XHq, with = 30 mm the height of the plate. 
The stiffness of the load cell is such that fracture occurs at fixed grips. A video movie of 
the plate is recorded at a typical 15 frames"^ rate and post-treated for tracking the crack 
tip position. Away from the sample edges, cracks run at a constant velocity V (see Fig.Q]). 
Further data processing is restricted to this region. V is computed by linear regression of 
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the tip position. Since reversible gels creep, the investigated range is restricted to velocities 
large enough for the energy released by stress relaxation of the gel to be negligible compared 
with that released by crack propagation. 

The energy released by unit area of crack advance is computed from the elastic energy of 
the uncracked plate ^(A) as^^ : Q = J^{X)/ (eL) whith e = 10 mm the thickness of the plate 
and L = 300 mm its initial length. This takes into account elastic non-linearities in the 
far-field region ahead of the crack tip. Each plate results in a single G{V) data point. We 
have also performed material-saving experiments in which the stretching ratio was increased 
at a constant rate A = 1.7 x 10^^ s~^. The resulting non-steady crack velocity along the 
crack path was computed from a piecewise linear fit. We have validated these Q{V) data by 
comparison with steady-state ones on an overlapping velocity range (Fig. jj^a)). 

Experiments with wetted crack tips are performed by injecting a drop of solvent of about 
250/il in the tip region while the crack runs. The solvent follows the tip and is prevented 
from flowing out both thanks to gravity (the crack running down vertically) and capillarity 
(the solvent wets the gel and forms a meniscus bridging the fracture gap). 
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FIG. 1: Velocity of a stable crack propagating in the mid-plane of a long plate. Sample: 
5 wt % gelatin gel in water, 10 mm thick, 300 mm long, stretched in the transverse direction by a 
factor A = 1.3. The crack was initiated by a single knife-cut notch, 20 mm long. Edge effects extend 
over about 2ho with ho = 30 mm the width of the unstretched sample. The remaining central part 
behaves as a homogeneously stretched plate, hence the observed steady crack-propagation. 




5 10 15 20 25 30 

V (mm/s) 



9 



120 
10D 
80 
60 
40 

20: 





1^ (wt %) 


11 (inPa.s) 









1 


3.5 


1.1x10' 


20 


1.8 


4.1 


1.5x10' 


30, 


2.5 


4.2 


1.6x10° 


60 


11 


S.,2 


2.0x1 0' 



□a' 



10 20 30 40 
riV (iiJ/m^) 



50 



m 



60 



FIG. 2: Influence of solvent viscosity on the fracture energy G{V). a, Gels with a constant 
gelatin content (5 wt %) in water /glycerol solvents with glycerol concentration ^ = wt % (circles), 
20 wt % (triangles), 30 wt % (squares), 60 wt % (diamonds). Filled symbols correspond to 
stationary cracks, open symbols to cracks accelerated in response to a steady increase of A. Qo = 
2.5 ± 0.5 J.m~^ is the common linearly extrapolated toughness, b, Plot of Q vs. r] V. Data points 
have been randomly decimated for clarity. Note the weak systematic growth of dQ/d{rfV) with 
glycerol content ^ elastic modulus /x (see Table). 
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FIG. 3: Effect of tip homo-wetting. G{V) curves for a 5 wt% gelatin gel in pure water : "dry" 
cracks opening in ambient air (upper data) and "wet" cracks with a drop of pure water soaking the 
tip. Each pair of points for a given Q corresponds to a single sample. Error bars show standard 
deviations of velocity along the track. At Q too low for dry cracks to propagate, wet ones can still 
run. Linear fits are shown. The wet data appear merely translated towards lower energies. The 
extrapolated fracture energy for wet tips is Q^^^ = 0.6 it 0.15 J.m"^. 
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FIG. 4: Effect of tip hetero-wetting by a less viscous solvent. The tip of a crack propagating 
through a gel of gelatin in = 30% glycerol/water is wetted with pure water (squares). The full 
curve is a guide for the eye. At low velocities, the data fall on the curve for a gel of gelatin and pure 
water wetted by pure water (full line). At higher velocities (^i < V2), the GiV) curve crosses 
over and approaches the one for the glycerolled gel wetted by the same solvent (closed circles and 
dash line). The data for a = 30% glycerol/water gel fractured in ambient air (open circles) are 
shown for comparison. 
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